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ABSTRACT 

This paper reports the measurement and prediction of three dimen- 
sional flow field in an axial flow inducer operating at a flow coefficient 
of 0.065 with air as the test medium. The experimental investigations 
included measurements of the blade static pressure and blade limiting 
streamline angle, and measurement of the three components of mean velocity, 
turbulence intensities and turbulence stresses at locations inside the 
inducer blade passage utilizing a rotating threads ensor hotwire probe. 
Analytical investigations were conducted to predict the three-dimensional 
inviscid flow and to approximately predict the three-dimensional viscid 
flow by incorporating the dominant viscous terms into the exact equations . 
Radial velocities are found to be of the same order as axial velocities 
and total relative velocity distributions indicate a substantial velocity 
deficiency near the tip at mid-passage® High turbulence intensities and 
turbulence stresses are concentrated within this core region. Evidence 
of boundary layer Interactions, blade blockage effects, radially inward 
flows , annulus wall effects and backflows are all found to exist within 
the long, narrow passages of the inducer, emphasizing the complex nature 
of inducer flow which makes accurate prediction of the flow behavior 


extremely difficult. 
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NOMENCLATURE 


C 

C f 

F 

h 

LE ,TE 


P 


Local blade chord 
Skin friction coefficient 
Viscous forces 
Static head 

Leading and trailing edges respectively 
Static pressure 


Pressure surface 

Total relative velocity (in experimental results non- 
dimensionalized with respect to U^) . 

Non-dimensionalized radius (= r/r ) 

Rotating cylindrical coordinate system (Fig. 1) 

Reynolds number (Q^ r/v) 

Local radius 

Coordinates (See Fig. 1) 

Suction surface 

Circumferential distance from suction surface non- 
dimensionalized by local blade spacing (distance 
between pressure and suction surfaces) . 

U»V,W Mean relative velocity components in the r, 0, Z 

directions respectively (in experimental results non- 

dimensionalized with respect to U, ) Fig. 1 

b 

U fe Blade tip speed (= r t ^0 

u,v,w Fluctuating relative velocity components in the r, 0, 

Z directions respectively (in experimental results, 
non-dimensionalized with respect to Q-) 

X Distance along the chord measured from leading edge 

Vg Absolute tangential velocity 

a Blade limiting streamline angle (Fig. 1) 

Flow angle (Fig. 1) 


Pr 

R 

r ,6 5 Z 
R 

e 

r 

s,n,r 

Su 
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NOMENCLATURE (continued) 


P 


3 

q 


Fluid density 
Shear stress 

2 

Static head coefficient (= 2g h/U^ ) 
Angular velocity of inducer 


Superscripts 


Time- averaged or passage-averaged quantity 


Subscripts 

t 

r,6>z 
s ,n,r 


Refers to values at the inducer tip 
Components in r, 8, z directions (Fig. 1) 
Components in s, n, r directions (Fig. 1) 



INTRODUCTION 


The Investigations reported in an earlier paper (1) is continued to 
obtain an accurate knowledge of the three dimensional mean flow and turbu- 
lence characteristics inside the passage of an axial flow inducer. The 
measurements, carried out with a rotating hot wire probe, and analysis, 
based on the exact solution of the equations of motion, has given a great 
insight into the flow characteristics, hitherto unknown, in this as well 
as other types of turbomachinery. 

THEORETICAL ANALYSIS 


Exact Inviscid Analysis ; 

Cooper and Bosch (2) have developed a method of obtaining the exact 
inviscid solution of the inducer flow field. Their three-dimensional 
analysis employs an iterative numerical procedure to solve the equations 
of motion expressed in finite-difference form. Cooper and Bosch program 
was subsequently modified by Poncet and Lakshminarayana (3) to include the 
trailing edge condition and by Gorton (4) to improve the convergence time. 
It has also been extended to include a viscid solution capability based 
on empirically determined blade skin friction coefficients. 

The nonlinear partial differential equations governing the flow in a 
rotating cylindrical coordinate system r, 0, z (Fig. 1) are: 
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In the Cooper-Bosch method, the above equations are rearranged to 
give residuals which are reduced to zero by a relaxation procedure. The 
total residual (RT) of one relaxation cycle is calculated by 


IMAX 

JMAX 

KMAX 

— 
RT - / > 

z 

7. + (R2 > 2 + (R3 ) 2 + (R4) 2 ] i)jjk (5) 

i=l 

j-1 
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where Rl, R2, R3 and R4 are the residuals calculated for the three momentum 
equations (Equations 1 to 3) and the continuity equation, and TMAX, JMAX 
and KMAX are the number of grid stations in the radial tangential and 
axial directions which are used in the numerical analysis. From this, the 
total RMS (root mean square) residual is defined as 


RMS = 


(RT) 

4 • (IMAX) • (JMAX) • (KMAX) 


( 6 ) 


and is thus a measure of the degree of convergence between the iterated 
solution and the exact solution. The boundary condition to be satisfied 
on the hub, annulus walls and the blade surfaces is ♦ n - 0 , where n 
is the direction normal to the channel boundaries and is the total rela- 
tive velocity. 

In applying Cooper-Bosch program to the Penn State inducer, the flow 
is assumed to be incompressible, and a grid of 7 x 7 x 26 is chosen to 
represent the blade passage. The flow geometry is shown in Fig. 1. 
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The first of the 26 axial stations corresponds to the upstream through flow 
boundary where the initial conditions are applied. The last four axial 
stations correspond to the downstream flow-through boundary, and extend 
to about one-fifth of the chord length downstream of the trailing edge. 

With Q 9 n = 0 to be satisfied on these stagnation stream surfaces, the 
set of boundary conditions for the problem is complete. 

For increased efficiency, the program written by Cooper and Bosch (2) 
has been compiled under a Fortran IV H level optimization procedure which 
reduces the time required for repetitive calculations. The next approach 
for the speedier solution of the governing flow equations is the optimiza- 
tion of the input parameters of velocity and pressure which would allow 
faster convergence to the three-dimensional solution. Cooper and Bosch 
(2) have derived an approximate solution to be used as an initial input to 
the exact program. This method derives the blade— to-blade average quanti- 
ties using axisymmetric equations, then uses these quantities in a blade- 
to-blade solution of an integrated form of the scalar momentum equation in 
the tangential direction. The flow parameters derived by this method were 
used in Ref. 3. However, an alternate method of developing the initial 
Input parameters is to estimate the velocities and static pressures inside 
the passage using the two dimensional solution. The method adopted here is 
to use the Douglas-Neutnann program described in Ref. 5. The two dimensional 
analysis is modified to include the effect of converging walls using the 
analysis of Ref. 6. 

Using the input parameters of velocity and pressure derived from the 
preceeding analysis results in a lower total RMS (root mean square) residual 
than with the previous method of initializing the input variables. As an 
example, the final RMS residual for the inviscid results of Ref. 3 was 
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0.12450 after 68 relaxation cycles, whereas a similar value is obtained 
using the present analysis in 10 relaxation cycles. This amounts to a 
considerable saving in computer time. Twenty-five iteration cycles has 
reduced the RMS residual to 0.10579, indicating that a faster convergence 

to the solution should be possible. 

In a further attempt to decrease the convergence time, the exit flow 
angle was allowed to change depending upon the tangential and axial veloci- 
ties calculated at the inducer trailing edge. Since the exact downstream 
boundaries are not known in this type of problem, it was hoped that by 
allowing the downstream boundaries to adjust themselves and thereby 
unload the blade trailing edge, a more exact definition of the downstream 
streamlines would result in lower RMS residuals. Cooper and Bosch suggest 
a similar technique as a means of reducing RMS residuals in their recom- 
mendations for future work. 

gifted the extension of the stagnation stream surfaces downstream have 
been constructed to be uniformly periodic with a spacing of 2 tt/N (N being 
the number of blades) , the values of velocity and pressure at the downr 
stream tangential channel boundaries should be equal. This condition is 
applied at the blade trailing edge after each iteration cycle. If the 
pressure and suction surface parameters differ with each other at the 
trailing edge grid point, the average value is used in the residual calcu- 
lations. If the axial and/or tangential velocities at the trailing edge 

diverge significantly from the design values during the iteration process, 

-1 V 

then the flow exit angle, defined by 3 — tan ^ at the trailing edge, is 
recalculated and is used to redefine the downstream stagnation stream 
surfaces. This method also has the advantage of automatically forcing the 
Kutta-Joukowski condition for the blade pressure distribution to be 
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satisfied. Changes made to the original Cooper-Bosch program including a 
listing of the computer program is given in Ref. 4. 

Cooper and Bosch program was run for the three bladed configuration 
(Fig. 1) without the modification' as well as with the modification forcing 
the pressure to be equal on both the pressure and suction surface at the 
trailing edge. The iterations were carried out for 24 cycles. The blade 
static pressures at the hub and tip radial locations are shown in Fig. 2. 
The modification close the pressure diagram perfectly. The predictions 
agree well with the experimental data. The hub wall static pressures were 
measured midway between the two blade surfaces and the annulus wall static 
pressures were derived from the static pressure taps located at various 
axial locations along the annulus wall. Hence the latter measurements 
are time averaged values. 


Approximate Viscid Analysis : 

In addition to the attempts to improve the convergence of the exact 
inviscid solution, a method of incorporating viscid effects into the 
governing equations of motion has also been investigated. 

In the momentum equations used by Cooper and Bosch (Equation 1 to3) , 
the following expressions for F^, and (the exact program variables 
for viscous loss terms) can be given as : 
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Molecular viscosity terms have been neglected in these equations. 

Since the stagger angle is very large, these viscous terms can be 
approximated by retaining the diminant terms as well as neglecting the 
normal shear stresses, resulting in: 
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The distribution of shear stress is assumed to be linear across the 
flow passage from pressure surface to suction surface (with zero at the 
midpassage) . The values of wall shear stresses are assumed to be known 
from previous experimentation. Skin friction coefficient for a four- 
blade flat plate helical channel is given in Ref. 7. The results, summarized 
in Fig. 3 are considered to be valid for the three-blade inducer under con- 
sideration. Interpolation of the curves in Fig. 3 for a given blade surface 

grid location under consideration gives a value of wall shear stress = 

1 —2 — 

* r p Q for the appropriate Reynolds number R = CL. * r/v of the flow 

I / K “Is. 

at that point, where is the average relative velocity across the flow 
passage as derived by the Cooper-Bosch relaxation procedure. The resultant 
turbulence stress at each grid point are thus calculated. The components 

i 

along each coordinate direction is assumed to vary according to the ratio 
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of the corresponding local velocity component to total relative velocity. 

The derivatives are calculated by finite difference method according to 
equation 10. 

An additional requirement placed on the viscid analysis is to satisfy' 
the viscid boundary condition which requires that all components of velocity 
are zero at the blade surface. The changes to the original Cooper and 
Bosch exact program necessitated by the inclusion of the viscous loss terms 
are made in subroutines "MAIN 1 ' , ”DLOSS*' and "RESXD" of the Cooper-Bosch 
program and are given in Ref. 4. It should be emphasized here that this 
is a first attempt to solve the entire viscous equation for rotor passage 
and is a preliminary one in as much as the shear stress distributions are 
assumed. For more accurate analysis, one has to resort to turbulent energy 
equation with suitable turbulence modelling. 

The input variables and formats for the modified viscid analysis pro- 
gram are identical to the original Cooper-Bosch program, with the exception 
of including a set of curves to define blade skin friction coefficient 
(C f ) vs. Reynolds number (R ) for various reference tangential locations 
throughout the inducer channel. Preliminary running of the modified viscid 
program indicates an increase in computer time of approximately two to 
three times more than a corresponding invlscid analysis run. This increase 
in computation time is due to the calculation of the viscous loss terms 
F , Fq and F^ at each grid point location throughout the duration of one 
relaxation cycle, which may involve several thousand iterations of the 
flow parameters in order to reduce the RMS residual from the previous 
cycle. 

The Cooper-Bosch exact analysis program incorporating the modifications 
mentioned in this section was run for the three-bladed Penn State inducer 
geometry. Both inviscld and viscid cases were considered. The inviscid 
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program was run. for approximately 75 relaxation cycles, resulting in a 
total RMS residual of 0.077. The viscid program was run for approximately 
50 cycles and produced a total RMS residual of 0,200. Most of the contri- 
bution to the total RMS residual comes from the upstream and the leading 
edge stations. In the viscid program the residuals at the blade surfaces 
were large in addition to the above stations. At most other locations, the 
total residual (local) was verying from 10 to 10~ . The results of the 
inviscid and viscid analyses will be discussed and compared with the experi- 
mental results later. 

EXPERIMENTAL FACILITY AND PROGRAM 

The experimental Investigation was performed in a 0.915 m diameter 

axial flow inducer with three equally spaced blades. Design of the rotor 
and description of the facility is given in Ref. 8. The use of the three 
bladed inducer for the continued experimental investigation is a result of 
the conclusions reached in Ref. 1; namely it has appreciably better per- 
formance than a four bladed inducer tested at the same operating conditions. 
The inducer was operated at a flow coefficient of 0.065 at 450 RPM with air 
as the test medium. The Reynolds number based on tip radius was 7 x 10^. 

An extensive experimental investigation of the rotor blade static pres- 
sure distribution has been undertaken to help provide a check on the theo- 
retical analysis and useful information for future theoretical development. 

A scanivalve, three channel pressure transfer device, slip ring unit were 
used to transmit pressures from rotor blade to a stationary manometer. 

The schematic of this set up is shown in Fig. 4. The blade static pressures 
were measured along quasi two dimensional stream surfaces at five radii 
shown in Fig. 5. The blade static pressures were measured at ten chordwise 
locations each on the blade pressure and suction surfaces. The coordinates 

of the measuring stations are tabulated in Ref. 4. 

The measurement of the blade limiting streamline angle (a) (Fig. 1), 
yhich is the limiting position of the streamline as the blade surface is 
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approached, is a further attempt to define the flow phenomena within the 
blade boundary layer, including nature and magnitude of the radial flows 
and the direction of the wall shear stress. The Information gained will 
help establish the extent of three-dimensionality in the inducer flow and 
will be valuable in developing a viscid theory for the prediction of rotat- 
ing boundary layer characteristics. The blade measurement stations are 
identical to those used for the blade static pressure measurements. The 
method and equipment for measurement are essentially the same as that used 
in Ref. 9. 

A triple sensor hot wire probe (rotating with the inducer) is used 
inside the passage to measure the three components of mean velocity, turbu- 
lent intensities and stresses. The governing equations, method, technique, 
instrumentation used and an estimate of the error involved are described 

fully in Ref. 10. The resultant voltage measurements from the hot wire 

~2 

were converted to mean velocities U, V, W and turbulence quantities u , 

~2 ~~2 — — — 

v , w , uv t uw and vw from the hot wire calibration curves and the appli- 
cable hot wire equations (Ref. 10). Measurements were taken at two axial 
stations, corresponding to approximately 33% and 90% of the blade chord 
(Fig.’ 5). Various velocity measurements have been performed at these 
stations (Ref. 1) and thus a comparison of hot wire experimental results 
with these prior investigations are possible (Ref. 4). Six radial stations 
(R = 0.973,. 0.945, 0.890, 0.781, 0.671, 0.548) at axial station A* and five 
radial stations (R = 0.973, 0.945, 0.890, 0.781, 0.671) at axial station 
B were traversed at several tangential intervals within the blade passage 
in an attempt to get an accurate and detailed appraisal of the flow 


Axial stations referred to as A, B and C in this paper are same as stations 
1, 2 and 3 in Fig. 1 of Ref. 1. 
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velocities, turbulence intensities and stresses in these regions. It was 
not possible, however, to obtain measurements close to the blade surface 
due to limitations caused by the blade curvature. 

EXPERIMENTAL RESULTS AND COMPARISON WITH THEORY 
Blade Static Pressure : 

Experimental and theoretical results are plotted in Fig. 6 for the 
five radial passage locations defined previously. It should be reiterated 
that the measurement stations do not correspond to constant radii, since 
the annulus area Is continuously varying. The measurement stations are 
shown in Fig. 5 and follows approximately the streamlines calculated from 
one dimensional consideration. 

The inducer design characteristic of trailing edge loaded blade is 
apparent from the measured ^ distribution. The blade static pressures 
are negative only near the leading edge of the suction surface and this 
follows the design trend based on cavitation criteria. The measured values 
of are found to be considerably different from the design values. The 
decrepancy is found to be maximum at the tip radius as shown in Fig. 6e. 

The design values are the ones assumed for deriving the blade profile from 
Wislicenus's (Ref. 8) 'mean streamline method.' Even though the design 
incorporates the effect of change in axial velocity and the blade blockage, 
the use of correlations based on cascade results are inadequate. This 
serves as a caution for those who are contemplating using the 'mean stream- 
line method' In designing unconventional blading. The experimental results 
indicate that the three dimensional inviscid effects are large. The flow 
near the tip, especially on the suction surface, is dominated by boundary 
layer effects (Figs. 6d and 6e) . 
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As can be seen in Fig. 6, the inviscid predictions agree well with 
experiment especially from hub to mid-radius, where the viscid effects 
are small due to large migration of the blade boundary layer towards the 
tip (Ref. 9). The inviscid and viscid predictions are close at most of the 
radii. The greatest variation between the two predictions occur near the 
tip (Fig. 6e) , where the viscous effects are dominant. The suction surface 
pressure distributions in this region predicted from the viscid analysis 
shows better agreement with the experiment than the inviscid analysis. 

In Fig. 7a are shown the passage averaged static pressure measured 

(using a stationary probe at the exit) very near the trailing edge (Ref. 3) 

compared with the predictions. Near the hub, the inviscid theory is in 

better agreement and at other locations (except at the very tip) , the 

viscid theory shows better agreement with the . experiment . The radial 

variation of ^ are compared with theory in Fig, 7b for axial stations A 
s 

and B shown in Fig. 5. The agreement between the theory and the experiment 
is very good, except at the very tip at station B. 

It is clear that the viscous effect has no appreciable influence on 
the blade static pressure distribution from hub to mid-radius. But its 
effect on velocities and shear stresses are considerable as indicated later. 

The measured and design static pressures are considerably different 
as shown in Fig. 7b for station A. The descrepancy at station B is found 
to be much worse (Ref. 4). 

Limiting Streamline Angles : 

The variation of the limiting streamline angle (a) with blade chord 
for each of the five radial measurement stations is shown in Fig. 8. The 
pressure surface distribution of a at the tip (radial station 5) indicates 
negative values of a (and, thus, radially inward flow) from leading edge 
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to mid-chord position. This tends to indicate the presence of the annulus 
■wall boundary layer scraping effect which induces flow away from the tip. 

At all other radial stations, a increases continuously from leading edge 
to trailing edge. Near the hub trailing edge, a increases quite rapidly. 

The blade limiting streamline angles at both radial stations 1 and 2 appear 
to extrapolate beyond 90° at the trailing edge, which is an indication of 
the existence of backflow in this region. This is presumably brought about 
by large radially outward flow that exists in the wake immediately down- 
stream of the trailing edge. This has a tendency to decrease axial velocity 
near the hub and thus induce backflows. At most axial locations, a 
decreases continuously from hub to tip. In several instances, this 
decrease appears linear. 

The suction surface a distribution remains relatively constant at all 
radial stations up to approximately 60% chord from the leading edge, when 
a more pronounced increase is noticed. At all stations except the tip, 
this increase extends to approximately 85%, chord and then a decreases 
toward the trailing edge. This is possibly due to the blade blockage 
effect in this region. At the tip, a increases continuously and no decrease 
is noted. Again, as in the pressure surface distribution, a decreases con- 
tinuously from hub to tip at practically all axial locations, and at some 
locations the variation appears linear. 

In most instances, the magnitudes of a on the suction surface are lower 
than at the corresponding position on the pressure surface. The deviation 
between pressure and suction surface measurements decrease continuously 
with increasing radius from hub to tip. The magnitude of ot, which is an 
indication of the extent of radial flows , is much higher than the values 
of a single blade reported in Ref. 9. This indicates that the radial 
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velocity in the inducers are quite appreciable, especially near the blade 
surfaces. 

Velocity Profiles at Station A : 

The axial, radial, relative tangential and total relative velocity 
components derived from hot wire measurements are described in this sec- 
tion. 

Fig. 9a shows the tangential variation of the total relative velocity 

(Q ) across the inducer passage at several radii. A distinct velocity 

deficiency is noted at approximately 55% from the blade suction surface 

for all radial stations, but it is especially pronounced near, the tip. 

This is the interaction region, where the radial flows inside the pressure 

and suction surface boundary layers, when encountered by the annulus wall, 

tend to roll towards mid-passage, interact, and produce strong vortices 

and radially inward flow as explained in Ref. 1. A concentration of high 

turbulence level in this region is confirmed by measurements described later. 

The tangential (relative) velocity profile, plotted in Fig. 9b, shows the 

same trend as Q_ thus indicating the dominance of relative tangential 
K 

velocity in inducers. Earlier results taken from a pitot tube (Ref. 1 - 
Fig. 12) do not indicate any deficiency in velocity at mid-passage. In 
view of the highly turbulent nature of the flow in this region, the pitot 
tube measurements may be in error at this location. The values of and 
V are found to be nearly same as design values from hub to mid-radius and 
considerably different from mid-radius to tip. The relative velocities 
near the tip region are much lower than design values resulting in large 
stagnation (absolute) pressure rise in this region. This large pressure 
rise is not due to flow turning but are caused by complex mixing and viscous 


interactions. 
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From the velocity profiles of Fig. 9a, it is easy to discern the edge 
of the suction surface boundary layer at radii above R = .671. The suc- 
tion surface boundary layer appears to grow in thickness as the tip is 
approached, increasing to approximately 25% of the passage width. This 
observation is consistent with the previous discussion about tip boundary 
layer interaction. No evidence of the pressure surface boundary layer can 
be detected in Fig. 9a. This tends to indicate that the suction surface 
boundary layer is thicker than that of the pressure surface, although it 
should be remarked that no measurements were taken close to the blade sur- 
face. Since the blade element is not radial, the hot wire probe could not 
be located very close to the blade surface. 

Fig. 9c shows the axial velocity (W) variation across the pas- 
sage width. The axial velocity is fairly uniform across the channel up to 
R 1=1 0.781. Beyond R - 0.781, it shows a tendancy to increase towards 

pressure surface. The radial variation of the axial velocity shows the 

largest values occurring near the hub, decreasing consistently towards the 
tip. This tends to indicate the effect of blade blockage on the axial 
velocity distribution. It is noted from Fig. 9c that negative values of 
W occur at the tip location R = .973. The existence of negative axial 

velocities at the extreme tip location indicates the presence of the annulus 

wall boundary layer scraping effect and was similarly noted by limiting 
streamline angle measurements at this location reported earlier. 

In Fig. 9d is shown the tangential variation of radial velocity U. 

Large values of radial outward velocity are found near the suction surface 
at radii close to the hub and negative radial velocities are found near 
the tip (R^O.781.), appearing at approximately 45% of the passage width. 

This is consistant with the previous discussion on boundary layer 
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interaction in this region. Maximum radial outward velocities are likely 
to occur near the blade surface and this trend can be discerned from Fig. 

9d. Near the hub, the radial velocities are induced mainly by flared hub 
and blade blockage, since the boundary layer is thin in this region as 
explained in Ref. 9. From mid-radius to tip, the radial flows are induced 
mainly by viscous and rotational effects. The radial velocities are of the 
same order of magnitude as the axial velocity, thus confirming the extent 
of three dimensionality in the flow. 

The total relative velocity distribution obtained from the viscid 
analysis is shown in Fig. 10a. The imposition of the boundary condition 
which defines the relative velocity on the inducer blade surface as zero 
enables the viscid program to provide a crude approximation for the pres- 
sure and suction surface boundary layers. The magnitudes of are similar 
to those found from the inviscid analysis, except near the blade surfaces. 

A slight velocity deficiency is noted near the tip at approximately 50/S 
passage width. This agrees with the experimental results of Fig. 9a and 
indicates an area of high viscous loss. Both inviscid and viscid analyses 
are compared with hot wire data in Fig. 10b for R = 0.973 and Fig. 10a for 
R “ 0.548. Away from the blade surface, prediction of Q R from -both the 
analyses are same and agrees reasonably well with the experimental data. 
Near the blade surface, the viscid analysis should provide better agreement 
with data. But there is no data available in this region to check the 
accuracy of the viscid analysis. It is also evident that the viscid analy- 
sis tends to over predict the extent of boundary layer growth near the hub. 
This is due to approximations involved in the viscid analysis as explained 
later. 

The axial velocities predicted from viscid and inviscid analysis at 
R “ .973 and R * .548 are shown compared with experimental results in 
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Fig. lOd. The predictions are found to be very similar at all radii. The 
predictions are good at R — .548 and poor at R = .973, thus indicating the 
dominance of viscous effects near the tip. No comments can be made with 
regard to the accuracy' of viscid results, since the measurements close to 
the proximity of the wall are not available. 

The radial velocities predicted from theoretical analyses are found 
to be very small in this region and the agreement between the theory and 
the experiment is only qualitative. 

Velocity Profiles at Station B : 

The tangential variation of the total relative velocity is shown 
in Fig. 11 for the axial station B and at various radial locations. Again, 
as in station A, a region of distinct velocity deficiency near the mid- 
passage is noted from mid-radius to tip. The velocity defect has increased 
considerably (at most of the radii) from that measured at station A. In 
comparison with the results of station A (Fig. 9a), the position of the 
maximum velocity defect appears to have shifted towards the suction surface 
to approximately 40% of the passage width. The measured values of Q are 
much lower than the design values at this location, except near the hub, 
indicating the extent of three dimensional lnviscid and viscid effects. 

The tangential variation of the axial velocity W is plotted in Fig. 

11 . Overall magnitudes are higher than those measured at‘ station A (Fig. 
9c) due to converging annulus. Again, as in station A, the maximum values 
occur near the hub radius and decrease > continuously towards the tip. The 
decrease in axial velocity near the tip indicates the continuing presence 
of the annulus wall boundary layer scraping effect. The effect, however, 
is not as severe as at station A, where negative velocities were measured 
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(Fig. 9c). It is also evident that the blade boundary layer increase in 
thickness as the tip is approached. 

The variation of radial velocity U across the passage is shown in Fig. 

11 . The overall magnitudes appear larger than at station A. The three- 
dimensional flow effects will be greater at station B and therefore 
accounts for the greater deviation of the flow from the two-dimensional 
design values which has been observed at this location. Fig. 11 indicates 
negative radial velocities foT the radii near the tip at approximately 25% 
from the suction surface. The radially inward flow at this location agrees 
with the previous discussions on boundary layer interaction and flow mixing 
which result in the velocity deficiencies and flow loss experienced in 
this region. 

The viscid analysis prediction for the total relative velocity distri- 
bution is shown In Fig. 12a. It is a striking departure from the inviscid 
analysis distribution (Ref. 4) especially near the tip where the viscous - 
effects are known to be appreciable. The viscid analysis also provides 
crude approximations for the suction and pressure surface boundary layers. 

A large velocity deficiency near the tip is predicted at approximately 50% 
passage width and agrees favorably with the experimental Q R profiles plotted 
in Fig. 11. . It is apparent that the velocity deficiency noted near the 
tip at station A (Fig. 10) has grown considerably as the flow proceeded 
downstream to station B indicating an increase in size and intensity of the 
vIscoub loss region and substantiating the experimental results discussed 
previously. The blade boundary layer development predicted by the viscid 
analysis can be seen in Fig. 12a. The suction surface boundary layer 
appears comparable in magnitude to the pressure surface boundary layer, 
both of which increase in thickness from hub to tip. The velocity gradients 
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near the suction surface boundary layer is steeper indicating the effect 
of adverse pressure gradients on this surface. 

Experimental and theoretical comparisons of distribution at R = 
0.973 and 0.548 are given in Figs. 12b and 12c respectively. It appears 
that the viscid analysis distribution provides better agreement with 
experiment, especially at the tip. The magnitude as well as the trend is 
predicted reasonably well. 

Fig. 12d shows the axial velocity distribution predicted from both 
inviscid and viscid analyses compared with the experiment at R = 0.671 and 
0.973. The viscid analysis reveals the approximate profiles inside the 
suction and pressure surface boundary layers. At R = 0.671 the predictions 
from viscid analysis agree better with experiment. While, near the tip, 
the predictions are poor from both the theories, especially near the pres- 
sure surface. The decrease in axial velocity from hub to tip is predicted 
and this is in agreement with, the data. 

The magnitude of radial velocities predicted from both inviscid and 
viscid analyses are found to be consistently lower than the experimental 
values. 

Passage Averaged Velocities : 

The passage averaged values of measured total relative velocity (Q^) , 
axial velocity (WT) and radial velocity (U) at stations A, B and at the exit 
are plotted dhd compared with the corresponding predictions from inviscid 
and viscid analyses in Fig. 13. It is evident from Fig. 13a that the pre- 
diction of Q_ is poor at axial station A and reasonably good at station B. 

Iv 

There is considerable decrease in Q R towards the tip at station B and this 
is predicted more accurately by the viscid analysis. The descrepancy 
between the design values (not shown) and the experimental data is found 



19 


to be large, especially at station B, thus indicating the need to incor- 
porate the dominant three-dimensional invescid and viscid effects in the 
prediction of the flow in inducers. The descrepancy between design and 
measured values of Q 0 is as much as 50% (design Q = 0.88, measured Q = 

0.4 at R « 0.975) near the tip, a consequence of appreciable boundary layer 
growth and mixing in this region. 

The passage averaged axial velocities (W) at stations A and B are 
shown in Fig. 13b. The measured values show a rapid decrease in axial 
velocity towards the tip at station A. This trend is predicted by viscid 
analysis. The agreement between theory and experiment is very good at 
station B. The magnitude of W increases from station A to B, an effect of 
converging annulus. Again, as in station A, the distribution of W shows 
larger values near the hub indicating the continuing presence of blade 
blockage, which is large at this radial location. It is interesting to note 
that the opposite trend is observed in the measurements taken downstream 
of the trailing edge (Fig. 5 - station C) and plotted in Fig. 13d. As 
shown in Fig. 13d, the axial velocities are very low at the hub (back flow 
region) and very high near the tip. This implies that significant changes 
occur in the axial velocity profile as the flow leaves the inducer passage. 

These changes may be responsible, in part, for the back flow region shown 

/ 

in Fig. 13d. 

In Fig. 13c are shown the passage averaged radial velocities compared 
with viscid and inviscid analyses. At stations A and B, the measured radial 
velocities are of the same order of magnitude as the axial velocities, 
further confirming our conclusion that the flow in inducers in highly 
three-dimensional. The radial velocities are higher near the hub (mostly 
caused by blade blockage and flared hub) and decrease slightly towards 
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the tip. Furthermore, the average radial velocity increases continuously 
from leading to trailing edge. The predictions are qualitative at 
station B. One possible reason for the poor prediction at station B is 
that the real fluid effects are not accurately accounted for in our analy- 
ses. The turbulent stresses are found to be very high in this region, 
even away from the blade surfaces as discussed later, and this is not taken 
into account in the viscid analysis developed in this paper. 

In Fig. 13d, the passage averaged values of absolute- tangential velocity 
(V Q ) and axial velocity (W) at -the trailing edge location are compared 
With the experimental data at station C shown in Fig. 5. At the hub, the 
prediction of Vq is poor. Near mid-radius, the inviscid analysis seem to 
provide better agreement. The viscid analysis provides better qualitative 
agreement near the ti£. The steep rise in Vq towards the tip is predicted, 
qualitatively, by the viscid analysis. Axial velocities predicted by both 
viscid and inviscid analyses are almost similar (not shown). Here again 
the agreement between the theory and the experiment is good only in the 
middle third of the annulus. One cannot expect the axial velocities shown 
in Fig. 1.3d to be predicted accurately, unless the three-dimensionality and 
viscous effects in the wake are accurately simulated in the viscid analysis 
presented in this paper. 

Turbulence Intensities : 

The distributions of radial, tangential and axial intensities at station 
B, non-diraensionalized with respect to local total relative velocity, are 
shown in isocontour form in Fig. 14. The contours for all three intensity 
components are essentially the same, showing a "pocket" or "core" of high 
turbulence centered at approximately 40% passage width and R = .890. This 
coincides with the location of the maximum total relative velocity deficiency 
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noted in Fig. 11., The turbulence intensities are generally higher than those 
encountered in a stationary passage. The peak intensities occur in the 
mixing region near the tip* where the two boundary layers merge and 
generate considerable flow mixing. The flow energy dissipated during this 
process is responsible for the velocity deficiencies encountered near the 
tip region in Fig. 11. Another concentration of high turbulence is noted 
near the hub pressure surface and is an indication of the proximity to 
the pressure surface boundary layer. 

""2" 

The radial component of turbulence intensities / u are generally 
higher than those in axial and tangential directions. This observed 
behaviour is opposite to known behaviour in stationary passages, where 

~~2 *~2 ~Y 

the streamwise component is the dominant one, and v > w > u . It can be 

shown, using the turbulent energy equation in rotating coordinates , that 

the rotation (or the coriolis forces) has marked effect on turbulence 

characteristics. Anand (7) has shown that in compressors and pumps, where 

the direction of rotation is opposite to that of the relative flow, the 

terms due to rotation yield an extra production term in the transport equa- 
~2 

tion for u , resulting in increased turbulence intensities in the radial 

direction. Anand* s (7) experiment also confirms the trend observed here, 

~~2 ~2 ~~1 
namely, u > v > w . 

It should be remarked that all the measurements reported here were 
taken away from the blade surfaces. The maximum turbulence intensities 
and stresses are likely to occur near the blade surfaces. High turbulence 
intensities measured away from the blade surfaces and reported here reflect 
the extent of turbulent mixing even near the mid-passage. The classical 
assumption that the viscous and turbulence effects are confined to very 
thin regions near the blade surface is evidently inapplicable to inducers 
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and hence, a fully three-dimensional treatment is needed for the predic 
tion of inducer flows. 


Turbulence Stresses : 

The distribution of turbulence velocity correlations vw, uw and uv 

were measured at station B using the method described in Ref. 10. What is 

required from the point of view of practical application is the shear 

stress in the cyl lnderical plane, V (In the direction of resultant 

velocity fv+W^) and radial direction, T , as shown in Fig. 1, The 

measured correlations (vw, uw and uv) can be used to find T and T as 

’ s r 

follows 


t * -pvv = (w cos.B + vsin$)(vcos3 - wsin$) 

8 11 , 9 


= (v 2 - w 2 )sinf3cos$ + vwCcos^P - sin 2 f3) 


( 11 ) 


where , 


sin$ = 


V + v 


cos(3 = 


W + w 



r (V+v) 2 + (W+w) 2 


/ 




( 12 ) 


(V+v) 2 + (W+w) 2 


v g and are fluctuating velocities in the direction of velocity 
(« V^^-Tw 2 ) and normal directions respectively (Pig. 1). 

Substituting equations (12) in (11) and neglecting second order terms 
compared to mean velocities we get 


T ~ — — w ^ VW (v — w 2 ) + vw(W 2 — V 2 ) L 
8 V^+W z l j 


( 13 ) 


similarly the turbulence stress in the radial direction is given by 
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* 


-puv - -p u(vcos|3 - w sin$) 
n 



( 14 ) 


The values of T , and |t{ so derived are plotted in Fig. 15. Con- 
centration of turbulence stresses are similar to turbulence intensity con- 
tours plotted in Fig. 14. The maximum stresses occur in the mixing region 
near the tip (from R = 0.78 to 0.973) at approximately 40% passage width. 

The values of T follows the same trend, qualitatively, as the velocity 

3 

profiles of plotted in Fig. 11 . Near the suction surface, away from 

the wall, the velocity gradients are negative resulting in negative value 

of stresses. Maximum stresses occur at R = 0.89 and near 30% passage width, 

this location corresponds to maximum velocity gradients (Fig. 11 ). The 

2 

turbulence stress (|t|/pQ^ ) measured near the wall of flat plate at zero 

- -3 

incidence is of the order of 1.5 x 10 . It is evident that the stresses 

measured (Fig. 15a) in a rotating inducer passage from mid-radius to tip, 
away from the wall, is of the same order of magnitude and is in fact 2 to 
4 times higher than this in some locations. This is one of the significant 
results of this investigation. More accurate estimate of the stresses 
within the channel is essential for the accurate prediction of the flow 
using viscid theory. 

The radial component of stresses are by far the most dominant ones as 
shown in Fig. 15b. As mentioned earlier, the stresses shown are away from 
the blade surface and their values near the blade surfaces are likely to 
be much higher. The radial stresses follow the same trend as the radial 
velocity gradients (note the correspondence between velocity gradients and 
T at R “ 0.973 and 0.945 at 20% passage point). One can prove qualitatively 
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(Ref, 7) , using Reynolds stress equation in rotating coordinates, that the 
effect of rotation is to increase turbulent stresses in the radial direction 
and suppress it in the streamwise direction. It is also evident that the 
total shear stress vector is not parallel to the resultant velocity vector. 
This is contrary to what is normally assumed in the analysis of three- 
dimensional turbulent boundary layer. 

f ~~2 2 

The resultant stress (x) (== /x +x ) are plotted in Fig. 15c. The 

S 17 

2 2 

ratio (|x|/pq ), where q is the total turbulence energy is found to vary 
from 0.04 to 0.27, with values averaging 2 in the region of higher stresses. 
In some boundary layer computations, the turbulence modelling is based 
on the assumption that the ratio of (stress/intensity) is constant, which 
is true only for isotropic turbulence. That this is not true is evident 
from this investigation. 


DISCUSSION AND CONCLUSIONS 

It should be remarked here that the viscid analysis carried out in this 
paper is a first attempt ever to solve for the real fluid effects in turbo- 
machinery passage using the entire equations of motion. It is a logical 
extension of Cooper and Bosch program, which, we think, was a breakthrough 
in the three-dimensional inviscid analysis of turbomachinery flows . 

Further improvements and refinements of the viscid analysis presented in 
this paper are needed before the three-dimensional flow field can be pre- 
dicted accurately. It is evident from the data reported, especially the 
shear stress measurements, that the shear stresses are non zero even ,at 
mid-passage and the resultant stress vector is not necessarily parallel to 
the resultant velocity vector. If the analysis could be modified to incor- 
porate these effects, the prediction of velocity components, especially 
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the radial component (which is not predicted accurately in this paper) 
could be vastly improved. Some of the major conclusions of this paper 
are ; 

(1) The viscid analysis, even though approximate, provides better 
prediction of the three-dimensional flow field, especially near the trail- 
ing edge. More tangential grid stations would be needed to better define 
the shape of the three-dimensional boundary layer. 

(2) The limiting streamline angle measurements indicate large radial 
velocities near the blade surfaces increasing gradually from leading edge 
to trailing edge. Negative values of ct at the pressure surface .tip indi- 
cate radially inward flow due to annulus wall boundary layer scraping 
effect. 

(3) The measured blade static pressures agree with the theory, except 
near the suction surface of the blade tip. 

(4) Total relative velocity measurements indicate substantial velocity 
deficiency at the outer radius near mid-passage, indicating complex inter- 
action between blade boundary layer and annulus wall. This is responsible 
for high stagnation pressure (absolute) size measured at the exit of the 
inducer near the outer radius. 

(5) Inside the passage, the axial velocities are generally higher 
near the hub and lower at the tip. This trend is reversed as the flow pro- 
ceeds downstream of the trailing edge. A slight backflow is observed at 
the tip near the leading edge. 

(6) The measured radial velocities are found to be of the same order 
of magnitude as the axial velocity. A region of radially inward flow is 
found near the tip region at mid-passage. 

(7) Turbulence intensities are generally much higher than those 
encountered in stationary passages, especially near the tip. The radial 
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components turbulence intensity is found to be highest and this effect can 
be attributed to rotation, 

(8) The streamwise and radial components of turbulence stresses are 
found to be very high, especially in the tip region. The radial stresses 
are higher than the streamwise components and this effect can be attributed 


to rotation. 
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